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I.  INTRODUCTION 


,  ,  of  a  spinning  projectile  with  a  liquid  payload  is  sometimes 

briefer  than  planned.  The  liquid  develops  an  Intrlcatepattern  of  oscilla¬ 
tions;  If  the  frequency  of  one  of  these  oscillations  Is  close  to  one  of  the 
yawing  frequencies  of  the  shell,  the  yaw  can  grow  unsuitably  large.  The 
projectile  becomes  unstable,  its  flight  performance  deteriorates  rapidly  and 
It  falls  to  achieve  Its  Intended  range. 

Much  work  has  gone  into  attempts  to  understand  and  describe  mathemati¬ 
cally  what  goes  on  between  liquid  payload  and  projectile.  One  of  the  workers 
In  this  field,  C.  H.  Murphy,  recently  derived1  a  set  of  equations  deflnlnq— 
under  certain  conditions— the  transverse  moment  exerted  by  the  fully  spun-up 
liquid  on  the  projectile.  I  have  used  these  equations  as  the  basis  of  an 
Interactive  computer  program  for  generating  liquid  moment  coefficients. 

Murphy  presented  his  equations  in  a  form  that  may  be  Intellectually 
satisfying  but  Is  several  steps  short  of  programming  suitability.  In  par¬ 
ticular,  certain  obvious  tasks  (taking  Indicated  partial  derivatives,  obtain¬ 
ing  closed-form  expressions  for  Integrals,  etc.)  were  never  undertaken. 

aware,?f  th?se  omissions,  and  It  was  his  hope  that  I  would  present 
the  missing  equations  In  a  report.  This  is  the  report.  K 

mnmont6™!? f i  J 1  I  V  how  the  computer  program  obtains  the  liquid 

moment  coefficients.  The  remaining  four  sections  discuss  four  alternative 
tasks  that  the  program  will  perform  upon  request. 


II.  THE  TRANSVERSE  LIQUID  MOMENT  COEFFICIENTS 


In  Reference  1,  Murphy  studies  the  response  of  the  fully  spun-up  liquid 
wagnltude^0^  6  COn^n®  mo^on  constant  frequency  and  exponentially  changing 

1  ■  *  0s*  (2.1) 

where 


s  *  (c  ♦  1)  x  (2.2) 

♦  *  ft  (2.3) 

l  •  constant  (positive)  spin 

* 

K  ■  complex  constant 


1.  C.  .¥.  Murphy,  "Angular  Motion  of  a  Spinning  Projectile  With  a  Viscous 
Liquid  Payload ,"  Ballistic  Research  Laboratoinj ,  Aberdeen  Proving  Ground , 
Maryland,  BRL  Memorandum  Report  ARBRL-MR-06194,  August  1982.  (AD  A11867&) 
(See  also  Journal  of  Guidance .  Control ,  and  Dummies.  Vol.  6.  Julu- 
August  1936,  pp.  230-296.)  - -  J 


♦  \  v 


A  .VVVVV.VV 


and  where  X  Is  the  complex  yaw  in  an  aeroballlstic  non-rolling  coordinate 
system. 

The  parameters  In  Murphy's  complex  variable  s  have  the  following 
Interpretation: 


t  *  (coning  rate)/$,  the  nondlmensionalized  coning  frequency* 


c  ■  (yaw  damping  rate)/ (coning  rate),  the  nondlmensionalized 
yaw  damping  (e  <  0)  or  yaw  growth  (c  >  0)  rate.  (Thus  t  -  0 
is  the  Important  special  case  of  constant-amplitude  coning 
motion.) 


Murphy  considers  the  liquid  payload  to  be  confined  to  a  right  circular 
cylinder  of  diameter  2a  and  height  2c.  The  axis  of  this  cylinder  Is  colllnear 
with  the  projectile's  principal  axis  (the  spin  axis);  the  center  of  the  cylln- 

*  l00???- wa  <H,stance  h  forward  of  the  projectile's  center  of  mass.  The 
offset  h  will  be  taken  as  zero  in  the  body  of  this  report;  nonzero  h  is  rele¬ 
gated  to  Appendix  E. 

The  cylinder  may  be  fully  or  only  partially  filled  with  liquid.  In 
either  case,  the  liquid  Is  assumed  to  be  fully  spun-up.  For  the  partlally- 
fined  cylinder,  It  Is  assumed  that  the  liquid  has  formed  a  cylinder  of  diame¬ 
ter  2a  with  an  air  core  of  diameter  2b.  (Murphy  also  considers  the  case  of  a 
fully-filled  cylinder  with  a  central  rod.  That  Is,  the  free  surface  r  «  b  Is 
replaced  by  an  inner  lateral  surface  r  «  d,  where  2d  is  the  rod  diameter. 

This  central  rod  case  Is  discussed  in  Appendix  F.) 

Murphy  introduces  two  real  liquid  moment  coefficients,  C,cM  and  C, 
defined  by  the  moment  equation  LIM 


T  S-M  * 


(*.♦! 


where 


M;  y  *  t  N, ~  ■  the  transverse  liquid  moment 
in  the  aeroballistlc  system 


*Othe r  notation  appeal's  occasionally  in  the  literature:.  Some  authors  have 
used  the  symbil  x  to  denote  the  complex  frequency  that  would  be  written 
here  as  -is. 
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■  the  mass  of  the  liquid  In  a 
fully-filled  cylinder 

3  liquid  density 

ClM  3  ^SM  *  1  ^IM*  (2.6) 


^he  ^LSM  term  (2.4)  Is  the  liquid  side  moment,  a  moment  causing 

rotation  out  of  the  transverse  plane.  The  term  In  Eq.  (2.4)  Is  the 
liquid  In-plane  moment,  a  moment  causing  rotation  the  transverse  plane. 


The  bulk  of  Reference  1  Is  devoted  to  deriving  an  expression  for  Cj^, 

This  expression  Involves  earth-fixed  cylindrical  coordinates  x.  r.  e  and  fo 
complex  perturbation  variables  deflnea  In  Eqs.  (3.10-3.13)  of  Reference  1: 


ps  -  nondimenslonal  pressure  perturbation, 

us»  vs»  #s  *  nondimenslonal  velocity  perturbation 
components  In  the  x,  r,  o  directions, 
respectively.* 


Each  of  these  variables  is  a  function  of  x  and  r  and  each  -  In  Murphy's  linear 
theory  -  is  linear  in  K. 

Murphy  makes  the  assumption  that  each  of  the  four  variables  can  be 
expressed  as  the  sum  of  an  inviscld  and  a  viscous  term: 


Ps  (r.*)  *  Ps1  (r.x)  ♦  PS¥  (r,x) 
and  similarly  for  the  velocity  perturbations. 


*The  association  of  (u,  v,  u)  tcith  the  ordering  (x,  r,  6)  seems  logically 
sound  and  refloats  Murphy 's  adherence  to  a  standard  reference  on  letter 
symbols 2  and  to  the  notation  he  has  used  in  over  SO  publications  on  sym¬ 
metric  missite  dynmies.  Hooever,  other  tori  tors  on  liquid-filled  shell 
associate  (u,  v,  >j)  tot  in  (r,  0,  x)  {or  in  their  notation  (r,  0,  s))i  still 
others  associate  (u,  \>,  o)  oith  rectangular  coordinates .  This  nonuni  f oust ty 
is  only  mildly  amusing. 

“American  Standard  letter  Symbols  for  Aeronautical  Sciences , "  ASA  ¥10. ?- 
1954,  published  by  The  American  Society  of  Mechanical  Engineers,  October 
1954. 


U% 


Making  a  few  additional  assumptions  and  adhering  rigidly  to  linear 
theory,  Murphy  derives  the  following  relations: 

“  (CLM}  pi  *  ^CLM)pe  *  (CLMJvt  +  (CLMJve 


where 


(C,J 


IM'  "lac 
pi 


Q 

j  |  K'1  x  ps(a,x)  -  x2/a  J  dx 


«w 


pe  a 


~  ^  |x_1  r2  ps1  (r,c)  -  cr3/a2 


dr 


(C,M)  -  (2K  Re  ex)'1  t  [x^.la^]  dx 

vt  J-c  r  •  a 


(C^)  -  (K  Re  ar) 


1  fa  f  3<"sv  -  1  VS»> 


U 


3X 


rdr 


x  *  c 


(2.8) 

(2.9) 

(2.10) 

(2.11) 

(2.12) 


and  where  Re  *  $a2/v,  the  Reynolds  number. 

The  subscripts  on  have  the  following  meaning: 


Subscript 

P* 

pe 

vt 

ve 


Denotes  liquid  moment  due  to: 
pressure  on  the  cylinder’s  Uteral  wall 
pressure  on  the  end  walls 
vHcous  shear  on  the  Jateral  wall 
viscous  shear  on  the  end  walls 


The  next  three  sections  will  be  devoted  to  converting  Eqs.  (2.9  •  2.12) 
to  programmable  form. 

111.  THE  FUNCTIONS  Xk(x) 


Murphy  derives  expressions  and  conditions  for  the  perturbation  variables 
that  Involve  complex  functions  of  r:  Rk(r),  and  complex  functions  of  x:  Xk(x). 

In  this  section,  we  define  and  discuss  Xk(x)  and  some  related  constants. 
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8 


1M 


n 


vm 
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The  subscript  k  Is  the  so-called  axial  wave  number,  the  number  of  nodes 
In  the  liquid’s  axial  wave  pattern.  The  liquid  has  three  wave  numbers,  (k,  n, 
m):  non-negative  Integers  associated  with  the  axial,  radial,  and  azimuthal 
wave  patterns,  respectively.  Since  the  functions  In  Eqs.  (2.9-2.12)  are 
evaluated  for  m  *  1,  only  two  wave  numbers,  k  and  n,  concern  us  here. 


The  axial  wave  number  k  Is  either  zero  or  an  odd,  positive  Integer.  A 
term  from  the  k  *  0  mode  Is  always  present  In  to  satisfy  the  Invlscld  end- 

wall  boundary  condition.  However,  Murphy  gives  this  term  explicitly,  circum¬ 
venting  the  need  for  a  k  *  0  subscript  In  this  context.  An  additional  term 
for  the  zero  mode  is  required  when  the  center  of  the  cylinder  Is  offset 
axially  a  distance  h  from  the  center  of  mass  of  the  projectile;  this  case  has 
been  banished  to  Appendix  E.  Hence,  k  will  not  assume  the  value  0  In  this 
report: 

k  -  1,  3,  5,  ...  (3.1) 


Before  we  discuss  Xk(x),  we  Introduce  the  complex  ^  (s,  c/a.  Re), 

which  can  be  computed  as  the  solution  of  the  equation 


1  *  \  &c  tan  *k  *  0 


(3.2) 


where 


-(a/c) 


2  /(T 


LL.  f  J— 11  ♦  i(  =Y| 

♦  is)  L  /3  ♦  Is  Vl  -  Is/. 


(3.3) 


1  ♦  1 _ 

✓  2(1  ♦  is)  Re 


(3.4) 


We  are  not  concerned  here  with  the  derivation  of  these  and  subsequent 
equations;  we  are  merely  trying  to  show  what  oust  be  done  by  a  computer 
program  to  obtain  results.  One  of  the  first  tasks  for  our  program  is  the 
solving  of  Eq.  (3.2)  by  an  iterative  process  (the  Newton  method).  A  good 
first  estimate  (for  small  jj  |,  that  Is,  for  large  Re)  is  the  approximation 


Ml  -  «e)  * 

Murphy's  complex  functions  XjJx)  have  the  form: 


(3.5) 


*k(x)  •  sin  (i^  x/c)  . 


(3.6) 


Next, we  introduce  three  sets  of  constants:  v  V  v  all  dependent  on 
Xk(c).  The  first  two  are  needed  to  compute  ak  and  ak  is  needed  to  produce  the 
functions  R^r)  of  the  next  section.  We  have 


9 


(3.7) 


(3.8) 


where  j,  like  k,  takes  on  the  values  1,  3,  5,  ... 

The  complex  constants  are  obtained  as  the  solution  set  of  the  system 

E  V  *k  ■  V  i  ■  *•  3>  *•— *  (3-9) 

k  »  1 


where  the  value  of  N  is  somewhat  arbitrary.  In  practice,  N  In  our  program  Is 
set  at  29;  this  allows  k  to  take  on  fifteen  values.  In  general,  solving  sys¬ 
tem  (3.9)  Involves  Inverting  an  (H  *  l)/2  by  (N  ♦  l)/2  complex  matrix. 

It  should  be  noted  that  the  ak's  determined  In  this  manner  are  precisely 

the  least  squares  coefficients  when  x/c  is  approximated  by  a  truncated  series 
In  Xk(x): 


*  : 
c 


N 


k  -  1 


Xfe(x) 


In  particular 


(3.10) 


H 

£  \  Xk(c)  ■  I  (3.11) 

k  »  1 


Of  less  use  to  us,  but  still  Interesting,  is  the  fact  that 

N 

£  a^  5^  ■  2/3  .  (3.12) 

k  •  1 
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.S  *  %Ws  k  *  •*  %■-  ,v  »V 


1 


>iNiNaVA 


For  the  special  case  of  Infinite  Reynolds  number,  we  have: 


5a  *  6c  *  0 


* 

*k/2 

Xk(c)  ■ 

(-l)(k-1)/2 

(  0  If  j  *  k 

V 

(  1  If  j  »  k 

"  Jk  * 

[8/(fk)2]{-l) 

V  Re"1  =  0 


(3.13) 


IV.  THE  FUNCTIONS  Rk(r) 

Murphy's  complex  functions  Rk(r)  have  the  form 
y  r)  .  j\  Oj  <*k)  ♦  f„  Y,  (*k)] 

where 

Ak  *  Vr)  "  *k 

C(3  ♦  Is)  (1  -  1s)]l/2  ^ 

\  - - 

1  ♦  Is 

Is  the  Bessel  function  of  the  first  kind  of  order  1 
(we  will  only  be  concerned  with  1-0  and  1) 

Is  the  Bessel  function  of  the  second  kind  of  order  1 

and  where  Ek(s)  and  Fk(s)  are  complex  functions  defined  by  Eq.  (4.1).  This 
definition  sheds  little  light  on  how  Ek  and  Fk  can  be  evaluated.  By  deter- 
mining  conditions  that  Rk  must  satisfy  at  the  boundaries  r  ■  a  and  r  ■  b, 
Murphy  obtains  two  equations  in  the  two  unknowns  Ek  and  Fk: 


(4.1) 

(4.2) 

(4.3) 


wnere,  setunj  '  s,  U  -  o,  J *  nave* 

CU  !I«,!  V  (1  -  Is)  *  (1  *  Is)  Au2  la  T  u 

*  [1  *  is  -  (1  -  is)  ia]  A:  ,  ^ 

c12  ~  (fonu  identical  to  Cjj  with  J  replaced  by  Y  } 

C21  *  &  ('  ♦  *s)  *  (3  *  is)s2]  J,b  -  (1  .  is)  Akb 

c22  =  tfonn  identical  to  c21  with  0  replaced  by  Y} 

ci  -  2is  (1  ♦  is)  (3  ♦  is) 
c2  »  (b/a)s2  (1  »  is)2  (3  ♦  is) 

and  where 

*ka  ■  *k(»)  *  \  i/c 
J0i  ’  Jo  (Akj) 

Jla  •  J|  (As,) 

Y0.  '  ?0  <*u> 

*1*  ■  »1  (*k.) 
and  sinllarly  for  r  »  i». 

'  In  Refer*™*  77  *urpkt,  approrivtea  la  by  in  au  and  fki a  int*>- 

dwsaa  on  error  t*ae  U  negligible  asoept  at  very  Ioj  Reynold*  nw&ere. 


Solving  the  system  (4.4),  we  have 


r  _  C1  ‘  (c12/c22}  c2 

Ek  - ; -  (4.10) 


_  ^cll/c22^  c2  “  ^c21/c22^  C1 

k  - : -  (4.11) 


where 

Gk  3  C11  ’  <c12/c22>  C21  (4.12) 


The  special  situation  Gk(s)  =  0  will  be  discussed  In  Section  VII  on  eigen¬ 
values.  * 

Note  that  for  b  =  0  (that  Is,  for  the  usual  case  of  a  fully-filled 
cylinder),  c22  is  Infinite,  so  that  Eqs.  (4.10-4.12)  reduce  to 

Ek  a  cl/cll 

Fk  3  0  b  3  0  .  (4.13) 

Gk  *  CH 

Finally,  we  note  that 

r  Y  (r)  .  .k  Ak  [Ek  0„  (A,)  ♦  Fk  (Ak)j  -  R,  (4.14) 

r2  Rk‘  (r)  -  (1  -  V)  R,  -  r  R,'  (,.,5) 

where  the  primes  denote  differentiation  with  respect  to  r. 

V.  THE  TRANSVERSE  LIQUID  MOMENT  COEFFICIENT  EQUATIONS 

Expressions  for  the  perturbation  variables  in  terms  of  Xk(x)  and  Rk(r) 

are  given  in  Appendices  A,  B,  and  C.  These  expressions  are  not  themselves 
part  of  the  computer  program,  but  they  are  useful  in  the  derivation  of  the 
programmed  equations.  The  nature  of  this  derivation  is  outlined  In  Appendix 
0. 


> v*  v» ri-'.MVAi'-'.v,* 


.T  vv*"-;  vr?T'- 


In  the  body  of  this  report,  as  we  have  Indicated,  we  are  interested 
mainly  in  presenting  the  programmed  equations  for  Cj_$M  and  C|_ We  can  now 
do  so: 


(CLM* 


i  (c/a)2  f  is  (2  +  is) 


III' 

T  ^ 


!l+st  1 

2  a  5-1 


(5.1) 


i  j  S’  ['.•(!) 


,-ld 


(5.2) 


(C,  u) 


(c/a) ^  (1 


*  «»>  T,j 


(5.3) 


(c/a)  <5 /  1  ♦  is  \1/2f  /  b*V 

[T4  -(*>/•)  T4b  -  2  1.(lM,)(l.7)J  (5.4) 


where 


Ti  ■  £  Ma>  •  bk 


(5.5) 


T2  =  Z  a  Rk'(a)  *  bk 


(5.6) 


t3  *  Z  tRk(a)  '  a  Rjc‘(a)3  •  V2  •  Xk(c) 
T3b  •  Z  CRk(b)  -  b  Rk'(b)3  •  x^2  •  Xk(c) 


(5.7) 


(5.8) 


T5  " 


ZRk(«) 

•  Xk(c> 

£  Mb) 

•  yo 

l 

r  (i  ♦  is)  t,  ♦  i  r2 

1  -  is 

L  3  ♦  is 

(5.9) 

(5.10) 


(5.11) 


and  where  the  summations  are  for  k  ■  1,  3,  5»...N. 
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To  produce  tabulated  values  and  plots  of  and  Ck|M  versus  t,  the 
computer  program  carries  out  the  following  Instructions: 

1.  Accept  the  required  input: 

Re,  b/a,  c/a,  e  and  the  x  range  of  Interest. 


2.  Divide  the  x  range  into,  say,  several  hundred  equally-spaced  points 
Tl’  T2  *  *  *  * 


3.  Perform  steps  4-10  below  for  each  x^. 


1.  For  the  current  s  value,  compute  5.,  d_  [Eqs.  (3.4,  3.3)]  and 
Cj,  c2  [Eqs.  (4.7,  4,8)].  a  c 


5.  Perform  steps  5a  -  5c  below  for  k  »  1,  3,  b,...N  (*29): 


5a.  Iterate  on  Eq.  (3.2)  to  obtain 


5b.  Compute  Xk(c)  and  bk  [Eqs.  (3.6,  3.7)]. 


5c.  For  j  »  1,  3,  5,...N,  compute  bjk  [Eq.  (3.8)]. 


6.  Solve  the  system  (3.9)  for  all  the  ak's. 


7.  Perform  steps  7a  -  7g  below  for  k  »  1,  3,  5,...N  (-29) : 


7a.  Compute  \  [Eq.  (4.3)]. 


7b.  Compute  Aka,  Akb  [Eq.  (4.9)]. 


7c.  Compute  the  complex  Bessel  functions  at  r  -  a  and  b: 
JGa*  Jla»  Y0a»  Yla*  J0b»  Jlb*  Y0b»  Ylb* 


7d.  Compute  cjj,  Cj2,  c2j,  c22  [Eqs.  (4.5,  4.6)]. 


7e.  Compute  Ek,  and  Fk  [Eqs.  (4.12,  4.10,  4.11)]. 
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7f.  Compute  R^a)  and  Rk(b)  [Eq.  (4.1)]. 


7g.  Compute  aRk * (a)  and  bR^ * (b)  [Eq.  (4.14)]. 


8.  Form  the  sums  [Eqs.  (5.5  -  5.11)]. 

9.  Compute  (C,M)  [Eqs.  (5.1  -  5.4)]. 

pt*  pe,  vt,  ve 


10.  Compute,  print,  and  store  Cl$m  3  R  {C^}  and  C^j^  »  I  { Cy|}. 


11.  Plot  Cl5n  and  C^m  versus  x. 


!  am  occasionally  bemused  by  the  fact  that  the  entire  process  above  takes 
only  a  fraction  of  a  second  per  point  on  our  computer  (a  Digital  Equipment 
Corporation  VAX-11/780)  while  It  simultaneously  keeps  a  dozen  or  so  other 
terminal  users  happy. 


VI.  THE  LIQUID  ROLL  HOME NT  COEFFICIENT 

As  mentioned,  our  program  has  the  ability  to  perform  tasks  In  addition 
to,  or  Instead  of,  computing  C^  and  Cj_jM  versus  t.  The  simplest  of  these 

optional  chores  Is  the  computation  of  the  liquid  roll  moment  coefficient, 

Cj^,  versus  t.  Cy^  is  defined  by  the  moment  equation 

*vi  ■  t  i?i2clrm  f®-11 


where  is  the  axial  liquid  moment  (compare  Eq.  (2.4)). 

In  Reference  3,  Murphy  shows  that  for  his  linear  fluid  mechanics  assump¬ 
tions  C|_RM  is  vety  nearly  (or  exactly,  when  c  »  0)  the  negative  of  C^: 

^RM  *  'S-SM  *  'tc/2)  C1  *  <4/3>  (c/a>23  *  (6.2) 


Thus,  the  computation  of  Cy^  from  C^  Is  trivial  and  we  move  on  at  once  to 
more  Interesting  options. 


3.  C.  H.  Murphy,  "liquid  Payload  BoV  Moment  Induced  by  a  Spinning  and  Coning 
Projectile,"  U.S.  Ballistic  Beseanoh  laboratory ,  Aberdeen  Proving  Ground , 
Maryland,  PPL  Technical  Report  ARBRL-TF-02SB1,  September  198 3.  (AD  A15S682) 
(See  also  AIAA  Fap*r  85-2142,  August  1985.) 


VII.  EIGENVALUES 


For  a  given  Re,  a,  b,  c,  and  k,  the  complex  variable  Gk  of  Eq.  (4.12)  is 
a  function  of  s.  For  some  set  of  values  s  3  {s^},  G|<  will  *>e  zero: 


Gk  (skn>  =  0 


(7.1) 


These  values  of  s  are  by  definition  the  eigenvalues  of  the  system. 

Index  n  (the  radial  wave  number  mentioned  in  Section  III)  Is  a  positive 
Integer  and  k,  as  before,  takes  on  the  values  1,  3,  5,...  In  theory, there 
are  an  infinite  number  of  eigenvalues,  but  their  practical  Importance  (that 
is,  their  effect  on  a  projectile's  performance)  usually  decreases  as  k  and  n 
Increase. 

One  of  the  options  in  the  program  developed  from  Murphy's  equations  is 
the  determination  of  the  more  Important  of  these  eigenvalues.  This  Is  done  by 
solving  Eq.  (7.1)  by  the  Newton -iterative  method,  in  which 


^skn^i+l 


<skn'l  “ 


v<*> 


(7.2) 


5  *  <5k„>1 


where  the  prime  denotes  differentiation  with  respect  to  the  complex  variable 
s. 

The  value  of  k  must  be  specified;  the  value  of  n  is  implicit  in  the 
required  initial  estimate  (s^g.  For  a  given  c/a  and  k,  the  program  deter¬ 
mines  six  initial  estimates  (n»l,2, . 6)  by  interpolating  in  a  stored 

table.  This  table  contains  a  set  of  real  eigenvalues  (Re“*  »  0,  b  ■  0) 
versus  c/(ak)  for  each  of  the  six  n  values.  The  stored  range  of  Is  -1  to 

1;  a  portion  of  this  table,  0  to  0.6,  is  given  here  as  Table  1.  From  these 
six  initial  estimates  (the  value  "six"  is  arbitrary),  the  program  attempts  to 
find  the  first  six  eigenvalues: 


skl»  sk2»#**sk6 


by  iterating  on  Eq.  (7.2).  Success  here  depends  on  the  adequacy  of  the  esti¬ 
mates  from  the  stored  table.  These  estimates  necessarily  become  poorer  as  the 
given  Re  decreases  and/or  b  increases.  Hence,  it  is  possible  (but  rare,  in 
our  experience)  for  the  Iterative  process  to  converge  to  a  “wrong"  eigenvalue; 
that  Is,  to  an  eigenvalue  associated  with  some  other  n  than  the  one  indicated 
by  the  program.  (This  problem  could  be  avoided  -  or  at  least  passed  on  to  the 
user  -  by  a  program  change  that  would  require  the  user  to  Input  the  six  eigen¬ 
value  estimates.  So  far,  this  has  not  been  necessary.) 


v.v.sv.v..  .v... 


r.  «  1 


? 


3 


5 


6 


0 

.9949 

.01 

! . 0064 

.02 

i .0160 

.03 

1.0297 

.04 

1.0417 

.05 

l .0538 

.06 

1 .0661 

.07 

l .0786 

.OS 

1.0913 

.09 

1.1042 

.10 

1 .1173 

.11 

1.1306 

.12 

1  !  /./.  1 

A  •  A  *t-T  4. 

.13 

1.1579 

.14 

1.1720 

.15 

1.1862 

.16 

1.2008 

.17 

1.2156 

.18 

1.2307 

.19 

1.2460 

.20 

1.2617 

.21 

1  .2777 

.22 

1.2940 

.23 

1 .3107 

.24 

1.3277 

.25 

1.3450 

.26 

1.3628 

.27 

1.3809 

.28 

1.3995 

.29 

1.4185 

.30 

1.4379 

.32 

1.4783 

.34 

1.5206 

.36 

1.5652 

.38 

1.6123 

.40 

1.6621 

.42 

1.7149 

.44 

1.7710 

.46 

1 .8308 

.48 

1 .8947 

.50 

1.9632 

.52 

2.0369 

54 

2.1164 

56 

2.2027 

58 

2.2965 

60 

2.3991 

-4780  .3103 

■  .3144 

.‘*VOh  .3186 

.4968  .3228 

.5032  .3271 

•5098  .3315 

•5165  .3360 

.5233  .3405 

.5302  .3451 

.5372  .3498 

.5443  .3546 

.5516  .3595 

cent  i  /  /  r 

•7771  .JUH) 

.5666  .3695 

.5743  .3747 

.5822  .3800 

.5902  .3854 

.5983  .3908 

.6067  .3965 

.6152  .4022 

.6239  .4080 

.6327  .4140 

.6418  .4201 

.6511  .4264 

.6605  .4328 

.6702  .4394 

.6801  .4461 

.6903  .4529 

.7007  .4600 

.7113  .4672 

.7222  .4746 

.7449  .4900 

.7687  .5062 

.7939  .5234 

.8204  .5416 

.8486  .5608 

.8784  .5813 

.9102  .6031 

.9441  .6264 

.9803  .6514 

1.0192  .6782 

1.0610  .7070 

1.1061  .7382 

1.1549  .7721 

1.2081  .8090 

1.2661  .8493 
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/(ak) 

4 


.2291 

.1814 

.2322 

.1838 

.2353 

.1863 

.2384 

.1888 

.2416 

.1914 

.2449 

.1939 

.2482 

.1966 

.2516 

.1993 

.2550 

.2020 

.2585 

.2048 

.2621 

.2077 

.2657 

.2106 

.2695 

.2135 

.2732 

.2165 

.2771 

.2196 

.2810 

.2227 

.2851 

.2259 

.2892 

.2292 

.2934 

.2326 

.2977 

.2360 

.3020 

.2395 

.3065 

.2430 

.3111 

.2467 

.3158 

.2504 

.3206 

.2543 

.3255 

.2582 

.3305 

.2622 

.3357 

.2663 

.3410 

.2705 

.3464 

.2748 

.3520 

.2793 

.3635 

.2885 

.3757 

.2983 

.3887 

.3086 

.4023 

.3195 

.4168 

.3311 

.4323 

.3435 

.4488 

.3567 

.4664 

.3708 

.4853 

.3859 

.5056 

.4022 

.5275 

.4196 

.5511 

.4388 

.5769 

.4595 

.6049 

.4820 

.6356 

.5067 

.1501 

.1521 

.1542 

.1562 

.1583 

.1605 

.1627 

.1649 

.1672 

.1695 

.1719 

.1743 

.1767 

.1792 

.1818 

.1844 

.1870 

.1898 

.1925 

.1954 

.1963 

.2012 

.2043 

.2074 

.2105 

.2138 

.2171 

.2206 

.2241 

.2277 

.2313 

.2390 

.2471 

.2557 

.2648 

.2745 

.2848 

.2957 

.3075 

.3201 

.3337 

.3483 

.3642 

.3814 

.4002 

.4208 


The  computation  of  Gk ' (s)  deserves  a  comment  or  two.  It  Is  barely  possi¬ 
ble  to  derive  an  exact  expression  for  Gk';  in  an  atypical  fit  of  zeal,  I  have 

done  so.  The  difficulty  lies  in  the  fact  that  everything  but  the  temperature 
in  Newark  depends  on  s.  It  was  necessary  to  derive  expressions  for  6  6 

V*  V»  J0a'*  etc,»  etc*»  a  tortuous  chain  culminating  in  Gk'  as  a 
function  of  the  C^'s  and  their  derivatives.  I  won’t  Inflict  these  equations 

on  the  reader  for  two  reasons:  (1)  they  would  take  up  a  lot  of  space  and 
time;  (2)  they  aren't  really  necessary.  The  approximation 


v<*> 


Gk(s  +  as/2)  -  Gk(s  -  as/2) 

— 


(7.3) 


should  be  more  than  adequate  if  the  increment  as  is  chosen  with  a  little 
care.  This  approximation  eliminates  the  myriad  opportunities  to  err  that 
arise  in  deriving  and  coding  the  exact  expression. 


VIII.  THE  EFFECT  OF  THE  LIQUID  MOMENT  ON  THE  DAMPING 

In  Section  V,  we  showed  how  our  program  determines  Cl$m  (t)  and  CLjM  (t) 

for  a  fixed  e.  But  e  and  t  are  not  really  Independent;  they  are  related  by 
the  yaw  equation.  In  Reference  1,  Murphy  derives  the  following  relationship 
(in  slightly  different  notation): 


where 


e(  x) 


BI  CLSM  (c,t) 


2  b2  [1  +  Bj  CLIM  (e.x)]  -  1 


♦  e. 


B1  s  mLa  /:x 


(8.1) 


b2  -  i  -  <«sg> 


-1 


Ix  »  axial  moment  of  Inertia  of  the  projectile 


sg  »  gyroscopic  stability  factor 


and  where  the  aerodynamic  damping  term  is  -  for  the  purposes  of  our  program 
-  a  specified  constant. 


.vn-KFi-*  •»  *  vmbiw.  v>.v>  V*  V>V*.V>  VW* 


One  of  the  program  options,  then,  1$  the  determination  of  e(x)  from  Eq, 
(8.1)  for  a  specified  Bj,  82,  and  range  of  t.  Note  that  e  appears  implic¬ 
itly  on  the  right-hand  side  of  Eq.  (8.1);  thus  an  Indirect,  iterative  method 
of  solving  (8.1)  for  e  is  needed.  For  a  fixed  t,  an  e  value  Is  assumed;  C» 5^ 

( e, t)  and  Cj_ jm  (e,t)  are  then  computed  as  In  Section  V  and  a  new  e  obtained 

from  (8.1).  The  new  e  value  replaces  the  old  one  and  the  process  is  repeated 
until  (In  a  well-ordered  universe)  it  converges  on  the  proper  e(x).  This 
approach  Is  carried  out  for  a  set  of  t  values  over  the  specified  range. 

For  most  of  the  cases  we  have  considered  so  far,  (sg)"*  has  been  taken  as 
zero;  that  is,  B2  3  1.  The  parameter  B^,  however,  is  not  necessarily  a  con¬ 
stant.  The  user  has  a  choice:  constant  Ix  (and  hence  constant  Bj)  or  con¬ 
stant  transverse  moment  of  inertia,  Iy,  In  the  latter  case,  the  program 
computes  the  variable  lx  from  the  relation 


!x  M 


mLa 


^IM 


(8.2) 


To  specify  Bj,  then,  the  user  inputs: 


1.  a  and  c  (to  allow  the  program  to  compute  the  constant  mass  mj_  from 
Eq.  (2.5)); 

2.  either  Ix  or  Iy  (whichever  Is  constant). 


IX.  THE  PRESSURE  COEFFICIENT 

Murphy's  complex  pressure  coefficient,  Cp  exp  (i$p),  Is  a  measure  of  the 

fluctuating  part  of  the  invlscid  pressure.  In  Reference  1,  he  derives  the 
following  expression: 

is>(l).(t)-Sl  <r-x>l  <9-» 


where  S^(r,x)  Is  a  summation  defined  In  Eq.  (Al)  of  Appendix  A.  The  final 
optional  task  our  program  can  perform  Is  to  compute  the  absolute  value  Cp  and 
the  argument  ^  from  Eq.  (9.1). 

There  are  various  ways  In  which  this  could  be  done.  For  example: 

1.  Fix  e,  r  and  x;  compute  Cp  and  $p  as  functions  of  i. 
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2.  Fix  e,  t  and  r;  compute  Cp  and  ^  as  functions  of  x. 

3.  Fix  e,  t  and  x;  compute  Cp  and  $p  as  functions  of  r. 

So  far  we  have  had  a  need  to  encode  only  a  special  case  of  1  above. 
Namely,  for  x  a  c,  e  =  0  and  any  specified  r/a,  the  program  will  compute 

CP  (t>  e'*p<T)  ■  -(f) li  <2  -  *>  (f)  +  si  <r-c>)  <9-2> 


X.  SUMMARY 

The  program  based  on  Murphy's  equations  has  so  far  been  used  mainly  to 
compute  the  liquid  side  moment  coefficient  and/or  the  eigenvalues  Sj^. 

For  a  specified  Re,  b/a,  c/a,  e  and  a  range  of  t  values,  Is  computed 
from  Eqs.  (5. 1-5.4)  and  Definition  (2.6).  (The  same  equations  yield  C^jM  but 
Interest  in  this  In-plane  moment  coefficient  Is  low.) 

For  a  specified  Re,  b/a,  c/a  and  k,  the  elgenfrequencies  skn  are  computed 
for  n  =  1  -  6  from  Eq.  (7.1)  and  Definition  (4.12). 
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Ws1 

*  1  *s1  • 

W4  [<*/o  *,  -  h] 

(A10) 

*Sl 

-  1  Vs1  ■  • 

■  [<•/<■>  Sj  ♦  S2  -  2  (1  + 

is)2  (x/c)]  . 

(All) 

For  the 

first  partlals  with  respect  to  r,  we  have: 

a  3  ps1 
a  ~1F~ 

*  -  (c/a)  K 

:  [s2  -  (l  +  is)2  (x/c)] 

(A12) 

3us1 

3  3  r 

■ 1  *[rr 

7s  *  ^  ♦  1s)  ] 

(A13) 

a  3  Vs1 

a  3r 

■ 1  (f)£t 

(a/r)2  S.  -  (a/r)  S  (a/c)2  S.  1 

3  ♦  Is  '  1  V  Is  J 

(A14) 

a  3W$1 

.  ;r'(a/r|2  si  *  <a/r>  h  2  h' 

3r 

w  L 

3  ♦  Is  (1  +  Is)2  J 

(A15) 

*  ’  (ws1 

*  1  V 

/'c\ ;p2  (a/r)2  sl  *  *  <t/r>  h 

(»/c)2  (1  .  1$)  Sj 

](A16) 

3r 

L  3  ♦  Is 

(1  ♦  is)2 

*  J  ("st 

-  '  *S1> 

(«/c)  (3  ♦  Is)  K  $5 

ar 

(1  ♦  Is)2 

(A17) 

For  the  first  partial*  with  respect  to  x,  we  have: 


•Air  ■(})<■  [«•«>*({)•  Sl 


(A19) 


c 


1  K  S5 
1  ♦  IS 


3  v$i  1  (c/a)  K 
9x~  *  1  -  Is 


[;....)!  »  a/r>  S3  *  <1  *  tS>  S<  ] 


(A20) 


3  wc.  (c/a)  K 


r  2  (1  +  is>  (®/r)  si  +  2  sa  l 

[(i.h)2 - mH - iJ 


(A21) 


3  (*s1  +  1  v$1)  (c/a)  K 


-  [(a/r)  S3  -  s4] 


(*22) 


3  (w$i  -  i  v$i)  (c/a)  K 


[2  U  *  is)2  -{a/r)  S3  -  s4] 


(A23) 


APPENDIX  B.  VISCOUS  PERTURBATION  VARIABLES  AT  THE  LATERAL  HALL,  r  *  a 

In  Reference  l,  Murphy  gives  the  following  lateral  wall  boundary 
conditions: 

PsvU.x)  =>  2  5a  w$v(a,x)  (Bi) 

usv(a»x)  3  -usi(a,x)  -  1  (1  +  Is)  K  (B2) 

*sv(a,x)  -  4a{-*s1U.x)  -  [a  ApJ  ♦  1  (1  a  Is)  (x/a)  k)  (83) 

'  r  *»  a  ’ 

wsv(a,x)  =  -  w$1(a,x)  ♦  (1  +  Is)  (x/a)  K  .  (B4) 

At  r  a  a  and  for  <5a  *  0,  differentiation  of  these  variables  with  respect 

to  r  Is  equivalent  to  division  by  s  : 

a 


r  a  a  d  1 
wsv 

3r 


r  =  a 


*  2  wsy(a,x) 


usv(a»x) 


r  =»  a 


r  »  a 


v$v(«»x) 


r  »  a 


*syU.x) 


(B5) 


(B6) 


(B7) 


(B8) 
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APPENOIX  C.  VISCOUS  PERTURBATION  VARIABLES  AT  THE  END  HALLS,  x  »  *c 


For  h  =  0,  Murphy1  gives  or  implies  the  following  end-wall  boundary 
conditions: 


Psv(r,c)  =  0 


us»<r-c)  *  -5c  c  ST 


x  ■  c 


vSv(r,c)  3  -v$i(r,c)  ♦  i  (1  ♦  is)  (c/a)  K 


wsv(r,c)  «  -wsi(r,c)  ♦  (1  ♦  is)  (c/a)  K 


and  for  the  other  end  wall: 


Psv(r,-c)  3  0 


%(r.-c)  *  ♦  usv(r»0 


vsv(r.-0  3  -  vsy(r,c) 


«Sy(r.-c)  *  -  w$y(r,c) 


The  following  partial  derivatives  are  given  or  itaplied  in  Reference  1: 


x  *  c 


»MV>OOt 

l»  •CAM* _ 


7  era-raw  ’PrjPJi 


ft 

I 


I 

I 


a 


§ 

i 


l'&\ 


X  ■  c 


[=^] 


X  «  c 


[■>] 


X  *  c 


■l^-l 


X  »  c 


■  (1/2)  (a  -  B)  wsi(p,c) 


-  (1/2)  (a  ♦  B)  v$1(rtc) 


♦  1  (1  ♦  Is)  (c/a)  K  a 


-  (1/2)  («♦  b)  «si(r,c) 

-  (1/2)  (a  -  B)  vsi(r,c) 
♦  (1  ♦  Is)  (c/a)  K  b 


T  3  Kv  ’  1  Ol 

[c - ^ - J  -  -  B  Cwsi(r,c)  -  1  vs1(r,c)] 


X  *  c 


♦  2  (l  ♦  Is)  (c/a)  K  B 


(C7) 


(C8) 


(C9) 


(CIO) 


h 

$ 

i-' . 


where 


1/2 


a  ■  (c/a)  [-1  (3  +  Is)  Re] 


B  •  (c/a)  [1  (1  -  Is)  Re]^2 


(Cll) 


(C12) 


and  where  the  square  roots  are  such  that  the  real  parts  of  a  and  a  art 
positive. 
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APPENOIX  0.  DERIVATION  OF  THE  TRANSVERSE  LIQUID  MOMENT  COEFFICIENT  EQUATIONS 

Our  final  expression  for  (C^)  ,  Eq.  (5.1)*  is  obtained  by  substituting 

pi 

in  Eq.  (2.9): 

for  psi(a,x),  using  Eq.  (A6); 

for  PsvU.x),  using  Eqs.  (Bl),  (B4),  and  (A9). 


The  resulting  integrand  involves  xS,(a,x)  and  xS?(a,x).  Integration  follows 

from  the  definition  of  b,  :  1  c 

k 


-2 


^  x  Sj (a. 


x)  dx  =  T 


1 


(01) 


^  x  S2(a, 


x)  dx  »  T2 


(D2) 


Our  final  expression  for  (C^)  ,  Eq.  (5.2),  is  obtained  by  substituting 


pe 


for  ps^(r,c)  in  Eq.  (2.10),  using  Eq.  (A6).  The  resulting  integrand  involves 


r2Sj(r,cl.  Integration  follows  from  the  properties  of  Bessel  functions: 


S,(r,c)  dr  =  (c/a)2  [T3  -  (b/a)  T3b]  . 


(03) 


Oar  final  expression  for  (CLM)  ,  Eq.  (5.3),  is  obtained  by  substituting 

vt 

in  Eq.  (2.11): 


for 


[•^]  ■ 


using  Eqs.  (B8),  (84),  and  (A9); 


r  *  a 


for  a 


[•>]  • 


-jr-  J  »  Eqs.  (86),  (B2),  and  (A7). 

r  «  a 


The  resulting  integrand  involves  S3(a,x).  Integration  follows  at  once  from 
the  definition  of  S3: 


31 


tt. 


v> 


.-1 


j  S3(a,x)  dx  •  2T4  , 


(04) 


Our  final  expression  for  (CLH)  ,  Eq.  (5.4),  Is  obtained  by  substituting 

for  the  partial  derivative  term  In  Eq.  (2.12),  using  Eqs.  (CIO)  and  (All). 
The  resulting  Integrand  Involves  the  combination 


aSjfr.c)  ♦  rS2(r,c)  »  a£[rRk(r)]'  Xk(c)  . 
Integration  follows  at  once: 


(05) 


.-2 


L  tasi(r- 


c)  ♦  r  S  (r,c)]  dr  -  T.  -(b/a)  T 


4b 


(06) 
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APPENDIX  E.  NONZERO  AXIAL  OFFSET 


Here. we  consider  the  case  where  the  center  of  mass  of  the  liquid-filled 
cylinder  is  located  a  distance  h  forward  of  the  projectile's  center  of  mass. 


The  functions  ( r)  were  defined  for  k  =  1,  3,  5,...  by  Eq.  (4.1).  For 
nonzero  h,  an  R0  must  also  be  considered: 


Ro<r>  ■  G)  [  Y  *  Y] 

where  E0,  FQ,  and  6Q  have  the  same  form  as  it.  Eqs.  (4.10-4.12) 


(El) 


Eo  = 


C1  ■  ^c12/c22^  c2 


(E2) 


Fo  = 


(C1 l/c22^  c2  "  ^C21/C22^  C1 


(E3) 


Go  =  C11  "  ^c12/c22)  C21 


(E4) 


The  definitions  of  Cj  and  c2  [Eqs.  (4. 7-4. 8)]  are  unchanged,  but  the  c^j's 
reduce  to: 


Cjj  *  3  +  is 


12 


(1  4  26a)  (1  -  is) 


c2l  *  (b/a)  (3  ♦  is)  s< 


c22  *  (b/a)’1  (1  -  is)  (2  ♦  4is  -  s2) 


(E5) 

(E6) 

(E7) 

(E8) 


Note  that  for  a  fully-filled  cylinder  (b  *  0),  Eqs.  (E2,  £3)  reduce  to 
E0  «2  is  (1  ♦  is)  \ 


)  b  »  0  .  (E9) 

In  Reference  1,  Murphy  shows  that  the  effect  of  h  on  the  (C^)  equation 


Fo  ‘  0 


is  to  add  a  small  term: 


P* 
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y 


(E10) 


where 


and  where 


(C.  J  •  [RHS  of  Eq.  (2.9)]  +  H  . 
L”  pi  P* 


*  -  —  - . 

pt  2a2t 


h  +  c 
'h  -  c 


Cp0<a)  dx 


cPo<a>  ■  [r.(d  -  t»..P  ^ar.]f  _  a 
■(I)  [Eo  *  F0  - 1s  <2  +  1s']  • 


Hence,  we  have 


H  =  -  )-lh/-al-  [E  +  F  -  Is  (2  +  Is)] 

pi  T  oo 


For  b  *  0.  this  reduces  to 


(Ell) 


(E12) 


(E13) 


u  .  is2  {h/a)2 
V  - T - 


and  for  b  *  t  *  0, 


Hpl'  -*t  0»/»)2 


[b  -  0] 


[b  »  c  *  0] 


(E13a) 


(E13b) 


Offset  h  has  no  effect  on  (C[^)pe  or  on  (Cyn)vt«  The  effect  on  the 
remaining  component  of  can,  according  to  Murphy,  be  written  as 


(cLM)ve  »  CRHS  of  Eq.  (2.12)]  ♦ 


(E14) 


where 


hj _  f*  [c  i<ail3!a>r  '  h  *  C  r  «r 

T  J  b  *-K  **  J  x  ■  h  -  c  ^ 


(E15) 


■  -  *  -  *  - ,  i  *  ^  «  'n.  vj  rK.vvwi.':7^  '>T>;*J<’ 


From  relations  given  by  Murphy,  we  can  simplify  the  integrand: 


3  (wsy  -  ivsy)  1  x  “  h  +  c  40  (h/a)  [EQ  -  21 s  (1  +  is)] 


x  =  h  -  c 


(E16) 


where  0  is  defined  in  Eq.  (C12)  and  where 


0  Re"1  =  (c/a)  (1  +  i) 


\  2  Re  I 


(£17) 


Substituting  (E16)  and  (E17)  in  (E15),  we  have 


Hve  =  (h/c)  ^C/a^ 


/  £  \  U  +  i)  [E0  -  2is  (1  +  is)] 

\  a2  /  t  [2  (1  -  is)  Re]1/2 


(E18) 


Thus,  the  effect  of  offset  h  on  the  liquid  moment  coefficients  can  be 
written  as: 


CLSM  =  (CL$M)  +RjHpi  +  Hve} 


(E19) 


CLIM  =  (CLIM*  +  1  (Hpi  +  Hvet 


(E20) 
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APPENDIX  F.  CENTRAL  ROD 


If  an  air  core  of  radius  b  is  replaced  by  a  central  rod  of  radius  d  the 
formulas  for  Ek  and  Fk  change.  System  (4.4)  becomes  ’ 


C11  Ek  +  c12  Fk  *  C1 


(FI) 


C11  Ek  +  c12  Fk  =  (d/a^  ci 


where 


C11  =  [(1-2  0(1  -  is)  -  (1  +  is)  A2  6  1  J 


kd  a  uld 


+  a  -  1s  +  (1  -  is>  sa]  Akd  J0d 


(F2) 


c12  =  Form  identical  to  c^  with  J  replaced  by  Y 
and  where* 

6  *  a  (d  +  a  O'1 


Akd  =  Ak  <d>  =  xk  d/C 


J0d  a  Jo  <Akd> 


r  (F3) 


Jld  =  J1  (Akd} 


and  similarly  for  YQd,  Yld#  The  rodded  Ek  and  Fk  values  are  the  solutions  of 
system  (FI),  the  rodded  determinant  Gk  Is  given  by 


(GkVod  =  C11  c12  "  ^11  C12 


(F4) 


and  the  rodded  eigenvalues  are  the  roots  of  the  equation 

Wrod  a  °* 


(F5) 


In  Reference  2,  Murphy  approximtea  Sfl  by  a  5  a/d  in  ou  and  ou.  Thie  ia 
oonaiatent  with  hia  uae  of  fia  for  5a  in  Eq .  (4.5)  but  ia  Uaa  valid.  The 
diatinotion  between  (d  +  a  6a)’J  and  d'1  oan  be  eignifioant  for  emit  rod 
radina  d. 
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The  rodded  value  of  CLM  win  depend,  of  course,  on  the  Ek * s  and  Fk  1  s 

computed  from  Eq.  (FI)  rather  than  Eq.  (4.4).  In  addition,  CLM  Is  affected  by 

explicit  changes  In  Eqs.  (5. 1-5.4).  For  the  two  end-wall  components,  the 
change  consists  of  replacing  b  with  d: 


where 


[(CLM)pe 


the  {CLM)pe  of  Eq.  (5.2)  with  b  replaced 
by  d  and  T3b  replaced  by  T3(j 


[<CLM>ve]  =  the  (Clm)V6  of  Eq.  (5.4)  with  b  replaced 
ro<i  by  d  and  replaced  by 

T3d  ■  £  £Rk<d>  -  dRk<d>}  •  ;k2  •  xk<c> 


T,d  *  E  Md>  •  Xk<c> 


For  the  two  lateral  components,  new  terms  appear  In  Eqs.  (5.1)  and  (5.3): 


i'v„l  ■  ^{('  -  S)  ^ 


-  7  CTt  -  (d/a)  Tld] 


(F10) 


where 


*  sa  tT5  *  T6d^ 


[<ClM>vtl  *  -  '-T-  j  T4  *  <d/»!2  T< 


*  ft  -  CT5fTM]j 


(Fll) 


Tid 3  L  Vd>  •  *>k 


( F12) 
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radius  of  the  liquid-payload  cylinder 

coefficients  In  the  least-squares  fit: 
x 


£  ak  Xk(x)  a  c 


[k  *  1,3,5, ...N] 


radius  of  the  air  core  In  a  partially  filled  cylinder 


-2 


C 


Xk(x)  x  dx 


[k  =  1,3,5,...] 


c"1  ^  Xj  (x)  Xk  (x)  dx 


[j,  k  =  1,3,5,...] 


half-height  of  the  liquid-payload  cylinder 

radius  of  a  central  rod  within  the  liquid-payload  cylinder 


axial  offset:  the  distance  between  the  centers  of  mass  of 
the  projectile  and  Its  liquid-payload  cylinder 


axial  wave  number  [1,3, 5,...] 

azimuthal  wave  number  (taken  as  1  In  this  report) 


2a  a^c  p^,  the  mass  of  the  liquid  In  a  fully  filled 
cylinder 


radial  wave  number  [1»2»3,...] 

nondlmensional  pressure  perturbation 


invlscld  component  of  ps 


viscous  component  of  ps 


radial  coordinate  In  an  earth-fixed  cylindrical  system 


(e  ♦  1)  t 
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LIST  OF  SYMBOLS  (Continued) 


gyroscopic  stability  factor 


eigenvalue  of  s  for  wave  numbers  k  and  n 


time 


axial,  radial  and  azimuthal  components  of  the 
nondlmenslonal  velocity  perturbation 


us1’vs1*ws1  1 "viscid  components  of  u$,  vs,  ws 


usv’vsv,wsv  v1scous  components  of  us,  vs,  ws 


axial  coordinate  In  an  earth-fixed  cylindrical  system 

A 

A^r/c,  the  argument  of  the  Bessel  functions 

"H.®  /** 


1  -  <4  .,) 


-1 


liquid  In-plane  moment  coefficient,  Eqs.  (2.4,  2.6) 

Clsm  +  1  Clih»  transverse  liquid  moment  coefficient,  Eq.  (2.4) 

that  part  of  Cy*  due  to  pressure  on  the  end  walls 
of  the  cylinder 


(CLM)pi  that  part  of  Cy<  due  to  pressure  on  the  lateral  wall 
of  the  cylinder 


that  part  of  Cy<  due  to  viscous  shear  on  the 
end  walls  of  the  cylinder 

that  part  of  Cy<  due  to  viscous  shear  on  the  lateral 
wall  of  the  cylinder 
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Cp  nondlmenslonal  magnitude  of  the  complex  pressure 
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s„(*) 


Ski’) 


Ix,Iy 
d0»  J1 


Eqs.  (4.10  -  4.11)  [It  •  1,3,5,...]  K 

determinant  of  the  system  that  determines  E0  and  F0, 

Eq.  (E4) 

determinant  of  the  system  (4.4)  that  determines  tk  and 
Fk;  the  roots  of  Gk(s)  *  0  are  the  eigenvalues  * 
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axial  and  transverse  moments  of  inertia  of  the  projectile 
Bessel  functions  of  the  first  kind  of  order  0  and  1 

complex  constant  In  the  definition  of  ?,  Eq.  (2.1);  Murphy’s 
equations  mere  obtained  by  linearizing  with  respect  to  K 

rectangular  components  of  the  liquid  moment  in  an  aeroballlstlc 
non-rolling  system 


•rtjtrary  upper  limit  on  k  for  computational  purposes 
(taken  as  29  in  our  program) 


v'i  ff)[¥  •  ’4\ 


Reynolds  number, 
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*k  CEk  J,  (\)  *  Fk  Y1  (/y]  [k  .  1,3,5, ...3 
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summations  defined  In  Appendix  A  [1  *  1  -  5] 

summations  needed  to  compute  and  C^, 

Eqs.  (5.5  -  5.11) 

sin  (^  x/c)  [k  »  1,3.5,...] 

Bessel  functions  of  the  second  kind  of  order  0  and  1 

(c/a)  [-  1  (3  ♦  Is)  Re]1/2 

(c/a)  [1  (1  -  is)  Re]1/2 

(1  ♦  1)  [2  (1  ♦  Is)  Re]‘1/2 

1  f  3  ♦  Is  1  -  Is! 

rcm ?r  [s-  *  <T~ J 

(yaw  damping  rate)/(con1ng  rate) 

aerodynamic  damping  constant,  Eq.  (8.1) 

azimuthal  coordinate  in  an  earth-fixed  cylindrical  system 

nondimensional  frequency  number,  determined  from  Eq.  (3.2) 
Ck  *  1,3,5,...] 

[(3  ♦  is)  (1  -  is)]l/2  ^ 
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K  e  *  the  complex  yaw  In  an  aeroballlstlc  non-rolling 
coordinate  system 

density  of  the  liquid 

coning  rate/i 

eigenvalue  of  t  for  wave  numbers  k  and  n 

it 

spin  (assumed  positive  and  constant) 

orientation  angle  of  the  complex  pressure  coefficient, 
Eq.  (9.1) 

complex  conjugate 

derivative  with  respect  to  whatever  independent 
variable  is  present 

value  at  r  »  b 

value  at  r  »  d 
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